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ABSTRACT 



The frequency-independent RMS temperature fluctuations determined 
from the COBE-DyVK two year sky maps are used to infer the parameter 
Qrms-PS} which characterizes the normahzation of power law models of 
primordial cosmological temperature anisotropy. In particular, a 'cross'-RMS 
statistic is used to determine Qrms-PS for a forced flt to a scale-invariant 
Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7° and 
10° RMS temperature derived from both the 53 and 90 GHz sky maps, we flnd 
Qrms-PS = 17.0^2!! /wK when the low quadrupole is included, and Qrms-PS = 
19.4^2!! /wK excluding the quadrupole. These results are consistent with the n = 
1 flts from more sensitive methods (e.g. power spectrum, correlation function). 
The effect of the low quadrupole derived from the COBE-DMK data on the 
inferred Qrms-PS normalization is investigated. A bias to lower Qrms-PS is 
found when the quadrupole is included. The higher normalization for a forced 
n = 1 flt is then favored by the cross-RMS technique. As initially pointed out 
in Wright et al. (1994a) and further discussed here, analytic formulae for the 
RMS sky temperature fluctuations will not provide the correct normalization 
amplitude. 



Subject headings: cosmic microwave background — cosmology: observations 
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1. INTRODUCTION 

Analysis of the first year of results from the COBE-DMK (Smoot et al. f992; 
Bennett et al. f992; Wright et al. f992; Kogut et al. f992) unambiguously demonstrated 
the existence of the long sought-after cosmological anisotropy in the Cosmic Microwave 
Background (CMB). The observed anisotropy is consistent with that predicted by models of 
structure formation with power law initial fluctuations of gaussian distributed amplitudes 
and random phases. The subsequent analysis of two years of data from the COBE-DMK 
(Bennett et al. f994) has confirmed and refined the initial results. 

In principle, the observed sky-RMS on a given angular scale is a convenient number 
to use in the normalization of a particular cosmological model. Previously, Wright et al. 
(1994a) used the sky- RMS from the first year COBE-DMK sky maps smoothed to an 
approximate FWHM of 10°, (Ts/;j^(IO°), to determine the effective normalization Qrms-PS 
for the scale-invariant Harrison-Zel'dovich power spectrum, P{k) oc A;", where n = 1 and 
k is the comoving wavenumber. It was also demonstrated that it is essential to account for 
both instrument specific details, such as the exact beam response function (rather than 
using a gaussian approximation, for example), and data analysis specific details, such as the 
subtraction of the best-fit monopole and dipole from the maps, which perturbs the inferred 
normalization from that derived from standard analytic formulae. An analysis based on an 
integral moment, such as the sky- RMS at one particular smoothing, does not have sufficient 
power to discriminate between different cosmological models. However, if the above details 
are taken into account then it can indeed provide a useful criterion for the normalization 
of a particular model. It should also be noted that using the observed sky-RMS values at 
a number of smoothing angles could be considered a poor-man's power spectrum analysis, 
and can be used to attempt to distinguish between models (see Smoot et al. 1994). 

In this paper, we update and extend the analysis of Wright et al. (1994a) using the 
two year COBE-DMK data to infer the normalization for an n = I power law model. 
A cross-RMS between two maps is defined and derived either including or excluding the 
quadrupole. The cross-RMS is determined from the maps with no additional smoothing 
(which we shall refer to as the 7° smoothing, since the central lobe of the DMR beam 
is approximately described by a 7° FWHM gaussian), and after an additional 7° FWHM 
gaussian smoothing to 10° effective FWHM. A likelihood analysis is performed for each 
smoothing both individually and jointly (although the joint analysis is not used with the 
intention of distinguishing between different power law spectral slopes). In particular, we 
find that the inferred Qrms-PS is reasonably independent of data selection, and is consistent 
with those values obtained from more powerful techniques, such as power spectrum 
estimates (Gorski 1994; Gorski et al. 1994; Wright et al. 1994b) and the correlation function 
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(Bennett et al. 1994), when restricted to n = 1 models. 



2. 



DATA SELECTION AND SKY SIGNAL ESTIMATION 



The COBE-DMK two year sky maps are used in the present analysis. To minimize the 
contribution from galactic emission, only those pixels for which the galactic latitude | b | > 
20° are used. Residual high latitude galactic emission and systematic error contributions 
are ignored (but are expected to be small, see Bennett et al. 1992, 1994). The best fit 
monopole and dipole and, where noted, quadrupole are also removed from the cut sky using 
a chi-squared technique and uniform weighting. The monopole and dipole are not physically 
relevant for the inference of the normalization parameter of a given anisotropy model, 
whilst the quadrupole is the multipole most contaminated by residual galactic emission and 
systematic errors in the maps. For determining the sky- RMS at 10°, the data surviving 
the galactic cut are smoothed by a 7° FWHM gaussian kernel with uniform weighting. 
GET_SKY_RMS, described in Wright et al. (1994a), uses a similar smoothing kernel, but 
with weighting by the number of observations per pixel. The cross-RMS, ®_rm5, between 
two maps a and b is then defined by 



where the sums are over all pixels i surviving the galactic cut, and Wi is the weight assigned 
to that pixel. Smoot et al. (1994) have discussed the effects of various weighting strategies 
on higher order statistics in some detail. The main result is that, although weighting of 
order N^^ (where is the number of observations in a given pixel i) will minimize the effects 
of noise on the sky- RMS, such weighting will increase the effects of cosmic variance in the 
ensemble of simulated skies used to infer Q^ms-PS- Unit weighting is preferred here, since it 
is nearly optimal in the sense of minimizing the spread in the simulated RMS values, and is 
more easily compared to analytic techniques. 

The ®RMS is a statistically unbiased estimator of the true cosmological RMS 
temperature anisotropy. If we write the observed temperature in a given pixel z, T^, as the 
sum of a cosmic term, t^, and a noise term, cTj-, then for two maps a and 6, 



P^'-^ i i i i 

where Npi^ is the total number of pixels, and we have taken the cosmic term in both maps 
to be equal. The noise terms will approach zero when averaged over a large number of noise 
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realizations. However, in any given realization of the (x)_rm5, such as the actual data, there 
will be non- vanishing noise contributions. These can obviously influence the outcome of 
the Qrms-PS inference, and this issue is addressed below. Note that, under the definition 
of eq.(2), the (E)rms is equivalent to the estimator of excess variance defined by Wright 
et al. (f994a) and applied to the first year COBE-DMK 53 GHz maps. Its virtue is that 
it is more easily computed and more readily identified with the elimination of any noise 
contribution. 

Table I shows the observed (E)rms values determined from the independent first and 
second year COBE-DMK sky maps at 31.5, 53 and 90 GHz. The results are consistent 
with the presence of sky signal in both years of data. The lack of any frequency dependence 
of the form in the (E)rms (expressed in thermodynamic temperature units) without 
any correction for galactic emissions is consistent with a cosmic origin for the signal. 
Nevertheless, the 31.5 GHz channels are the most likely to be contaminated by residual 
galactic emission and are also less sensitive than the 53 and 90 GHz channels. For the 
remainder of the analysis we only consider data at the two higher frequencies. Table 2 
summarizes the (E)rms values derived from possible combinations of the 53 and 90 GHz 
A and B channel maps. Note that whilst the (E)rms values are generally in excellent 
statistical agreement, the outcome of the Qrms-PS analysis can be sensitive to the particular 
combinations selected. In the likelihood analysis below, three combinations are employed. 
The first, 53(A+B)(x)(90A+90B), is the one used in the correlation function analysis of 
Bennett et al. (1994), the second, 53A(x)53B, is the equivalent combination to that used in 
Wright et al. (1994a), and the third, (53+90)A(x)(53+90)B, is one of the map combinations 
used in Wright et al. (1994b). 

3. LIKELIHOOD ANALYSIS 

In this Letter we restrict our attention to the inference of the Qrms-PS normalization 
for the n = 1 power spectrum model. The analytical form of the probability distribution 
of the COBE-DMK data analysis specific (E)rms statistic is unmanagable. A Monte Carlo 
approach was adopted to generate the (E)rms distributions for a grid of Qrms-PS values 
(with 2500 simulations used for each value of Q^ms-ps)- Each simulation generates maps 
of the sky temperature distribution by combining a realization of Harrison-Zel'dovich 
sky anisotropy filtered through the COBE-DMK beam (Wright et al. 1994a) with noise 
realizations based on appropriate values of the rms per observation and observation patterns 
of the specific DMR channels. The effects of noise correlations are negligible (Lineweaver 
et al. 1994) and, hence, are not included in these simulations. The Qrms-PS-dependent 
statistical means, variances and covariances of the 7° and 10° (E)rms are derived from 
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these Monte Carlo simulations, and used to construct the gaussian approximation to the 
probability distribution of the (E)rms- This, together with the measured (E)rms values, 
defines the likelihood function C{Qrms-Ps)- 

Figure 1 shows the likelihood functions obtained from the data at 7° and 10°, and from 
a joint analysis of the two smoothing angles, both including and excluding the quadrupole. 
Table 3 summarizes the maximum likelihood values of Qrms-PS and the 68% and 95% 
confidence level intervals. These values of the rms quadrupole normalization amplitude for 
an n = 1 power law model as determined by the (E)rms technique are completely consistent 
with several previous analyses based on the two year 53 and 90 GHz data. Gorski et al. 
(1994) use a power spectrum analysis and find Qrms-PS = 19.9 it 1.6 fiK and 20.4 it 
1.7 fiK including or excluding the quadrupole respectively. Wright et al. (1994b) use a 
quadratic power spectrum estimator excluding the quadrupole to determine Qrms-PS = 
19.8 lb 2.0 fiK. In Bennett et al. (1994) the cross-correlation function is employed to infer 
Qrms-PS = 18.2 lb 1.5 fiK when the quadrupole is included, and 18.6 ±1.6 fiK excluding 
the quadrupole. 

It should be recalled that the first year 53 GHz 10° sky-RMS, analyzed solely with the 
quadrupole included, rendered the best estimate of Qrms-PS = 17.1 ib 2.9 fiK (Wright et al. 
1994a). The corresponding result from the two-year data is 17.3^2^1 (see Table 3), which 
agrees well with the one year result (with appropriately lower errors). However, our analysis 
of the combined 7° and 10° (E)rms with quadrupole excluded implies a higher Qrms-PS 
normalization. Other more powerful techniques as applied to the COBE-DMK data for 
the inference of this cosmologically interesting parameter have also proven sensitive to the 
inclusion of the observed quadrupole, which has a relatively low value of 6 ib 3 fiK (Bennett 
et al. 1994). Whether this diminuitive quadrupole amplitude has specific cosmological 
implications, is a result of chance cancellation of residual galactic emission with some of 
the cosmic quadrupole emission, or is only a manifestation of cosmic variance remains open 
to speculation. However, its consequences to this analysis are of interest, and have been 
studied with the aid of miscellaneous Monte Carlo simulations. 

The particular Qrms-PS = 20 fiK model selected for scrutiny was suggested by the 
analysis of Gorski et al. (1994). This was the most consistent, either including or excluding 
the quadrupole, in the inference of Qrms-PS for an n = 1 spectrum. Figure 2 shows a scatter 
plot of the (E)RMS at 7° and 10° smoothing from 2500 simulations of noise-contaminated 
Harrison-Zel'dovich skies. The (E)rms data point when the quadrupole is included is ~ Icr 
deviant from the ensemble average of the simulations (Figure 2a), which is statistically 
satisfactory. When the quadrupole is excluded (Figure 2b), the data and simulations are 
in excellent agreement. So, are these to be considered mutually exclusive? A tenable 
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argument to the contrary arises from additional simulations in which the realization-specific 
quadrupole is restricted to be less than 9 fiK. Figure 3a demonstrates the consistency 
between the (E)rms data and the average over this ensemble. As seen in Figure 3b, this 
constraint does not affect the quadrupole-excluded case. This, then, provides a plausible 
explanation for the difference in Qrms-PS amplitudes inferred when either including, or 
excluding, the small observed quadrupole from the analysis. To quantify this, a sample 
drawn from the restricted quadrupole simulations was processed by the maximum likelihood 
machinery. A bias of 2 - 3 fiK was observed to lower Qrms-PS (depending on which 
smoothing angles were involved) when the quadrupole was included. The bias corrected 
maximum likelihood Qrms-PS normalizations are 20. f, f9.5 and f9.7 fiK for the 7°, fO° 
and joint analysis respectively. No bias was observed when the test data were analyzed 
excluding the quadrupole. Thus, we can conclude that if the bias introduced by the low 
observed quadrupole is accounted for, then the normalizations inferred with the quadrupole 
either included or excluded are consistent in the context of a forced n = f fit. A Qrms-PS 
value of f 9.4 fiK is the most appropriate for the (E)rms analysis of an n = f spectrum, after 
accounting for the bias due to the low quadrupole. 

A further important aspect of the analysis demonstrated in Figure 2 is the comparison 
with analytic calculations of the RMS. The analytic formula for the sky- variance, with a 
gaussian approximation to the beam filter function, specified for an n = f model, is 



where £ is the spherical harmonic order and erf, is the gaussian beam dispersion. Here, erf, 
= 3° corresponding to the approximate 7° FWHM COBE-DMK beam. The sum over 
£ is taken in the range [2,40] or [3,40], which is sufficient to determine the (E)rms either 
including or excluding the quadrupole respectively. This formula overestimates the strength 
of the fiuctuations, since the actual COBE-DMK beam filter function drops in amplitude 
in £-space more rapidly than the gaussian approximation (see Wright et al. f994a). Fven 
when the correct beam and pixelization weights are included, there remains a disagreement 
between the analytic calculations and the simulations that explicitly account for the 
monopole, dipole and, if required, quadrupole subtraction from the galaxy-cut data. Since 
these multipole estimates are made on incomplete sky coverage, there is some aliasing 
of higher order power into the fitted and removed low order amplitudes, thus the RMS 
fiuctuation amplitudes are suppressed. The combined beam filter and multipole subtraction 
effects are of order 5-f 0% when the quadrupole is included (in agreement with Wright et al. 
f994a), and f0-f5% in the no quadrupole case. 
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4. DISCUSSION 

Some aspects of the analysis related to a proper assessment of the results include: 
1) choice of input data values, 2) noise model uncertainties, 3) biases in the parameter 
inference and 4) relation to other (power law) cosmological models. 

We have only used three of the cross-combinations possible with the 53 and 90 GHz 
data in this analysis, as motivated by previous work. Inspection of the observed (E)rms 
values in Table 2, together with the errors on the inferred Qrms-PS values, suggests that all 
of the combinations should be consistent with Qrms-PS ~ 19 /uK with the possible exception 
of the 90A(x)90B 7° RMS. However, this is most likely just an anomaly due to noise: in 2500 
simulations generated with Qrms-PS = 19 /uK, this particular combination yielded a zero 
®_RM5 in 7.0% of the simulations including the quadrupole, and in 12.2% of those excluding 
the quadrupole. Further, although the 7° (E)rms prefers a Qrms-PS normalization of zero, it 
is still consistent with the 19.4 fiK normalization determined previously at the ~ 2cr level. 
The 10° and joint analyses of the 90 GHz data, both when the quadrupole is included or 
excluded, are consistent with this normalization to ~ Icr. An analysis including both the 
53 and 90 GHz data is to be preferred due to its higher sensitivity. 

The noise RMS per observation, (Tobs^ is known to an accuracy of ~ 1%. Simulations 
performed with (Tobs adjusted by such an amount demonstrate that the maximum likelihood 
values for Qrms-PS are shifted by ~ 0.1 //K, an insignificant amount compared to the error 
in the inferred quadrupole normalization. 

The particular observed (E)rms value from the 53 and 90 GHz data is noise contaminated 
(eq.2). The important issue is to determine if this results in a biased inference of Q^ms-PS- 
A sub-ensemble of the simulations was selected and used as test input data. No statistically 
significant bias was observed in the simulated sample-averaged estimates of Qrms-PS either 
including or excluding the quadrupole. 

This analysis has been specific to a forced n = 1 spectral fit. The normalization of 
other models of cosmological anisotropy should proceed either by a detailed reworking of 
the above, or by using more powerful techniques that are sensitive to both Qrms-PS and n, 
such as the power spectrum method (Gorski et al. 1994; Wright et al. 1994b) or the 2-point 
correlation function (Bennett et al. 1994). In fact, the analysis of Gorski et al. (1994) has 
rendered a useful power spectrum independent normalization, ag ~ 8 //K, for power law 
spectral models, whilst the two-point correlation function technique described in Bennett 
et al. (1994) finds ar ~ 9.5 fiK. 

An exact calculation for the power spectrum of CMB anisotropy that includes all the 
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potential, velocity and adiabatic effects on the last scattering surface in the context of a 
primordial n = 1 spectrum renders an effective spectral slope slightly steeper than n = 1 
over the angular scales probed by the COBE-DMK instrument. However, since such an 
exactly computed power spectrum is not sufficiently described by a simple power law, we 
have implemented an analysis as above using multipole coefficients, kindly provided by 
Radek Stompor, for a specific h = 0.5, ftb = 0.03, CDM model. The Qrms-PS normalization 
inferred from the joint 7° and fO° analysis of the 53 and 90 GHz data is fG.Tti^g /wK 
including the quadrupole, and f8.9lj;g //K when the quadrupole is excluded. As discussed 
previously, the latter higher value is preferred in the present analysis for the low quadrupole 
case at hand. 

In summary, we have used the (E)rms statistic derived from the COBE-DMK two 
year 53 and 90 GHz maps to infer a normalization Qrms-PS ~ 19 /uK for the amplitude 
of primordial inhomogeneity in the context of a Harrison-Zel'dovich n = 1 model. We 
stress that simple analytic models that do not include a correct description of the 
COBE-DMK beam or the monopole, dipole and, if appropriate, quadrupole subtraction 
will underestimate the amplitude of Qrms-PS when normalizing to the observed sky- RMS. 
The low observed quadrupole amplitude affords a reasonable explanation for the difference 
in inferred Qrms-PS amplitudes when either including or excluding the quadrupole. 

We gratefully acknowledge the efforts of those contributing to the COBE DMR. 
COBE is supported by the Office of Space Sciences of NASA Headquarters. We thank 
Radek Stompor for providing us with the CDM model CMB anisotropy power spectrum 
coefficients computed using his Boltzmann code. 
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Table 1: Observed (E)rms in thermodynamic temperature derived Irom the first and second 
year 31.5, 53 and 90 GHz COBE-DMK data. The errors are determined from a large number 
of noise simulations for a fixed CMB sky realization. The frequency independence of the 
data, expressed as ®_rm5('^) is demonstrated. 



1st year (/) 2nd year 
®i?M5 combination 


including Quadrupole 


excluding Quadrupole 


7° (/.K) 


10° (/iK) 


7° (/.K) 


10° (/iK) 


31(A+B) 
53(A+B) 
90(A+B) 


39.7t^?:g 
35.3+t? 


QS + 7.6 
OO.Z_7 2 

32.7+11 


34.8t^^:° 

25.1^13.9 


33.2t?-^ 
30.9il| 

32.2t3;5 


Frequency dependence: /3 


-0.45 ± 0.8 


-0.1 ± 0.2 


-0.3 ± 0.9 


0.1 ± 0.2 



Table 2: Observed (E)rms values derived from possible 53 and 90 GHz combinations. Those 
combinations used in the Qrms-PS analysis are denoted by 'f. 



®i?M5 combination 


including 


Quadrupole 


excluding 


Quadrupole 




7° ifiK) 


10° (/iK) 


7° ifiK) 


10° (/iK) 


53A(g)53Bt 




32.4+1-^ 


43.lt^:^ 


30.6tl:^ 


90A(g)90B 




25.7t|l 




24.7t|^ 


53A(g)90A 


28.4t?:^ 


30.9t^:^ 


97 7-I-9.2 
^ ' • ' -8.6 


30.3+^-^ 


53B(g)90B 


09 0-1-7.0 


00 q + 2-3 


30.2t?:2 


27.7i^:^ 


53A(g)90B 


45.2t^i 


32.4t2;2 




30.9i^;2 


53B(g)90A 


34.6t|^7 


01 7-1-2.9 

' -3.0 


00 7-1-10.8 

'J'J- ' -10.0 


30.8t|° 


53(A+B)(g)90A 


31.6t^:g 


31.3t^i 


30.8t6':6 


30.6t^;5 


53(A+B)(g)90B 


39.3tl? 


01 9-1-1.9 

'^-'-•^-2.0 




00 Q-i-2.0 

^».<J_2.0 


53A(g)90(A+B) 


37.7t^i 


31.7il| 




3o.6i^:° 


53B(g)90(A+B) 


33.5i^-? 


30.8t^:l 


31.9iL' 


29.3t^:l 


53(A+B)(g)90(A+B)t 


35.7t^-° 


3i.2tl:? 


UAtti 


29.9tl:^ 


(53+90)A(g)(53+90)Bt 


35.6tt2 


30.7tl:^ 


34.4tl^ 


29.4tl:^ 


(53A+90B)(g)(53B+90A) 


30.311^ 


29.9tl:^ 


28.8it^ 


28.5tl:^ 
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Table 3: Derived parameters from the likelihood (C) analysis, assuming an n = 1 spectrum. 



53(A+B)(g)90(A+B) 

Qrms-PS (/^K) 


including Quadrupole 


excluding Quadrupole 


7° 


10° 


joint 


7° 


10° 


joint 


Maximum C value 


18.0 


16.6 


17.0 


20.2 


19.0 


19.4 


68% c.l. interval 


[14.7, 21.2] 


[14.7, 19.1] 


[14.9, 19.5] 


[16.6, 23.5] 


[17.1, 21.3] 


[17.3, 21.7] 


95% c.l. interval 


[9.8, 24.5] 


[12.9, 22.2] 


[13.1, 22.9] 


[10.9, 26.5] 


[15.2, 24.1] 


[15.5, 24.6] 




53A(g)53B 

Qrms-PS (/^K) 


including Quadrupole 


excluding Quadrupole 


7° 


10° 


joint 


7° 


10° 


joint 


Maximum C value 


22.3 


17.3 


16.7 


25.0 


19.4 


19.1 


68% c.l. interval 


[18.7, 25.8] 


[15.2, 19.8] 


[14.7, 19.0] 


[21.3, 28.6] 


[17.4, 21.9] 


[17.0, 21.4] 


95% c.l. interval 


[13.8, 29.9] 


[13.3, 23.1] 


[13.0, 22.3] 


[15.5, 32.4] 


[15.4, 24.8] 


[15.1, 24.3] 




(53+90)A(g)(53+90)B 

Qrms-PS (/^K) 


including Quadrupole 


excluding Quadrupole 


7° 


10° 


joint 


7° 


10° 


joint 


Maximum C value 


18.1 


16.6 


17.0 


20.4 


19.0 


19.4 


68% c.l. interval 


[14.5, 21.4] 


[14.7, 19.1] 


[14.9, 19.6] 


[16.4, 23.8] 


[17.0, 21.3] 


[17.3, 21.7] 


95% c.l. interval 


[9.3, 24.7] 


[12.9, 22.3] 


[13.1, 23.0] 


[10.3, 26.8] 


[15.2, 24.2] 


[15.4, 24.7] 
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Fig. 1. — Likelihood curves lor Qrms-PS derived Irom the (E)rms assuming n = 1. Thin 
lines represent the analysis based on the 7° RMS, thicker lines Irom the 10° RMS, and 
thickest lines Irom the joint analysis. Solid curves include the quadrupole, the dashed 
curves are for the no quadrupole case. Top: 53(A+B)(x)90(A+B), center: 53A(x)53B, bottom: 
(53+90)A(g)(53+90)B. 
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Fig. 2. — Scatter plots for 2500 simulations of the 7° and 10° (E)rms for Qrms-PS = 20 fiK 
and n = 1 with the noise properties of the 53 and 90 GHz (A+B) maps. The cross represents 
the 53(A+B)(x)90(A+B) data, the filled circle is the mean from the simulated ensemble, the 
square is the analytic prediction using the COBE-DMK beam and pixelization smoothing, 
and the diamond is the analytic prediction using a gaussian approximation to the COBE- 
DMR beam and no pixelization smoothing. Also shown is the Icr density contour derived 
from the simulations. Top: quadrupole included, bottom: quadrupole excluded. 
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Fig. 3. — As Figure 2, but with the ensemble comprising 2500 realizations with a realization- 
specific quadrupole less than 9 fiK. The cross represents the 53(A-|-B)(x)90(A-|-B) data, and 
the filled circle is the mean from the simulated ensemble. Top: quadrupole included, bottom: 
quadrupole excluded. 



